KAIST-TH 01/01, SNUTP 01-002, 



riep-ph/0101026 



Axino as a sterile neutrino and R parity violation 



Kiwoon Choi a , Eung Jin Chun 6 and Kyuwan Hwang a 

a Korea Advanced Institute of Science and Technology, Taejon 305-701, Korea 
1 'Department of Physics, Seoul National University, Seoul 151-747, Korea 

Abstract 

We point out that axino can be a natural candidate for sterile neutrino which would accommodate 
the LSND data with atmospheric and solar neutrino oscillations. It is shown that the so called 
3 + 1 scheme can be easily realized when supersymmetry breaking is mediated by gauge interactions 
and also imparity is properly broken. Among the currently possible solutions to the solar neutrino 
problem, only the large angle MSW oscillation is allowed in this scheme. The weak scale value of 
the Higgs fi parameter and the required size of -R-parity violation can be understood by means of 
spontaneously broken Peccei-Quinn symmetry. 
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Current data from the atmospheric and solar neutrino experiments are beautifully ex- 
plained by oscillations among three active neutrino species ffl. Another data in favor of 
neutrino oscillation has been obtained in the LSND experiment [0 . Reconciliation of these 
experimental results requires three distinct mass-squared differences, implying the existence 
of a sterile neutrino u s 0. In the four-neutrino oscillation framework, there are two possible 
scenarios: the 2 + 2 scheme in which two pairs of close mass eigenstates are separated by the 
LSND mass gap ~ 1 eV and the 3 + 1 scheme in which one mass is isolated from the other 
three by the LSND mass gap. It has been claimed that the LSND results can be compatible 
with various short-baseline experiments only in the context of the 2+2 scheme M. However, 
according to the new LSND results , the allowed parameter regions are shifted to smaller 
mixing angle, accepting the 3 + 1 scheme |5|, 0. Although it can be realized in a rather lim- 
ited parameter space, the 3+1 scheme is attractive since the fourth (sterile) neutrino can be 
added without changing the most favorable picture that the atmospheric and solar neutrino 
data are explained by the predominant — > v T and v e — > u^, v T oscillations, respectively. 
In particular, the 3 + 1 scheme with the heaviest v s would be an interesting explanation of 
all existing neutrino data. On the other hand, it is rather difficult to find a well-motivated 
particle physics model which would yield the desired four-neutrino masses and mixing in a 
consistent way ||. 

In this paper, we show that the 3+1 scheme can be easily realized in supersymmetric 
model with £7(1) Peccei-Quinn (PQ) symmetry when supersymmetry (SUSY) breaking is 
mediated by gauge interactions. In this model, axino can be as light as 1 eV, so plays the role 
of sterile neutrino |J. A proper axino-neutrino mixing can be induced by i?-parity violating 
couplings which appear as a consequence of spontaneously broken U(1)pq. It turns out that 
only the large angle MSW solution to the solar neutrino problem is allowed in our model. 
The weak scale value of the Higgs /^-parameter and the required size of i?-parity violation 
can be understood by means of the Frogatt-Nielsen mechanism |1(| of spontaneously broken 

U(1)PQ. 

The model under consideration includes three sectors: the observable sector, the SUSY- 
breaking sector, and the PQ sector. The first contains the usual quarks, leptons, and two 
Higgs superfields, i.e. the MSSM superfields. The second contains a gauge-singlet Goldstino 
superfield X and the gauge-charged messenger superfields Y, Y c . Finally the third contains 
gauge-singlet superfields Sk which break U{1)pq by their vacuum expectation values (VEV), 
as well as gauge-charged superfields T, T c which have the Yukawa coupling with some of Sk- 
The superpotential of the model is given by 

W = hS 3 (SxS 2 - f PQ ) + kSiTT c + \XYY C + H/mssm + (1) 

where PUmssm involves the MSSM fields, and Wsb describes SUSY breaking dynamics en- 
forcing X develop a SUSY breaking VEV: XX = M x + 2 F X . This VEV generates soft 
masses of the MSSM fields [|ll|, m so ft = -^F x /M x , which are assumed to be of order the 



weak scale. One can easily arrange the symmetries of the model, U(1)pq and an additional 
discrete symmetry to get 

^mssm= u!j 'I :'■>'■-; + V^xQiD] + y^QiU? + j^SfHxHz 

-SlLiH 2 + '^y-SiL i L j El + -^-SiLiQjDl, (2) 
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where the Higgs, quark and lepton superfields are in obvious notations. Here M* denotes 
the fundamental scale of the model which will be taken to be the grand unification scale, 
Mqut ~ 10 16 GeV. Full details of the model will be presented elsewhere [|12[ . 

Let A = (0 + ia) + 9a + 9 2 F A denote the axion superfield where a, and a are the 
axion, saxion and axino, respectively. It is then convenient to parameterize Si and 5*2 as 
Si = Se A ^ PQ and £2 = Se~ A ^ PQ . Note that the VEV of S is uniquely determined to be 
(S) = fpQ which would correspond to the axion decay constant if the saxion is stabilized 
at e^l$ p Q Ri 1. After integrating out the SUSY-breaking sector as well as the heavy fields 
in the PQ sector, the low energy effective action includes the Kahler potential and the 
superpotential of A, 

K eS = f PQ {e {A+Ar ^ + e-^y^} + AK cS 
W cS = fi e 2A '^HiH 2 + ^ A ^^UE 2 

+e A /^(\ ljk L t L J E c k + KikLiQjDl), (3) 

where AK e n is A-dependent loop corrections involving the SUSY-breaking effects and 

im> = yofp Q /M*, A = v'Jl Q /Ml 

Kjk = li 3 kfp Q /M„ X' ijk = j' ijk f PQ /M*. (4) 

This shows that in our framework the weak-scale value of fio and also the smallness of im- 
parity violating couplings can be understood by means of the Frogatt-Nielsen mechanism of 
U(1)pq with fpg M* [0. Although not written explicitly, the coefficients of B- violating 
operators UfD^D^. can be easily arranged to be small enough to avoid too rapid proton 



decay, including the decays into axino or gravitino [|T^]. The best lower bound on fpq is 
from astrophysical arguments implying fpq > 10 9 GeV [H|. To accommodate the LSND 
data, we need the axino-neutrino mixing mass of order 0.1 eV. It turns out that this value 
is difficult to be obtained for f PQ > 10 10 GeV. We thus assume f PQ = 10 9 - 10 10 GeV with 
M* = Mqut for which /i takes an weak scale value (with appropriate value of yo) and all 
i?-parity violating couplings are appropriately suppressed. 

An important issue here is the saxion stabilization. One dominant contribution to the 
saxion effective potential comes from AK c g which is induced mainly by the threshold effects 
of T, T c having the A-dependent mass My = ^fpQ^ A ^ PQ ■ If Mt ^ AX, one finds 



N T M T Ml ( A 2 Z T Z T . 



>-eff 



167T 2 Z T Z T c \ M T M, 



In "~'~r (5) 



T 



where Nt is the number of superfields in T, Z? is the Kahler metric of T, and A is a cutoff 
scale which is of order M X - If M T > XX, AK cS would be enhanced by \XX/M T \ 2 , however 
then the resulting saxion effective potential can not stabilize (f) at the desired VEV [12] . We 
thus assume M T < XX. With Z t \qiqi ~ — m 2 oit , the above AK eS gives a negative-definite 
saxion potential 

vi 1} ^-^rn 2 soit \Kf PQ \ 2 e 2 ^. (6) 
There is another (positive-definite) potential from the A-dependent /^-parameter: 

Vi 2) *e^^\^\ 2 (\Hi\ 2 + \H 2 \ 2 ). (7) 
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If nfpq is of order few TeVs, <\> is stabilized by + vj 2 ^ at the desired value e^^ PQ ~ 1. 
This requires a rather small Yukawa coupling k ~ 10~ 6 which may be a consequence of some 
flavor symmetry. For e^^ PQ ~ 1, the saxion and axino masses are estimated to be 

m J = (10 - 10 2 keV) 2 + AmJ, 

m~ a = (Kr 4 -l(T 2 eV)+Am a , (8) 

where the numbers represent the gauge-mediated contributions for f P q = 10 9 — 10 10 GeV, 
and Am^ and Am 5 are the supergravity-mediated contributions which are generically of 



order m 3 / 2 jnj. In gauge-mediated SUSY breaking models jnj, the precise value of m 3 / 2 
depends on the details of SUSY breaking sector. However most of models give m 3 / 2 > 1 eV, 
implying m„ is dominated by supergravity contribution. In this paper, we assume 

m a ~ Am a ~ 1 eV (9) 

which would allow the axino to be a sterile neutrino for the LSND data. We note that 
although it is generically of order m 3 / 2 , Am 5 can be significantly smaller than m 3 / 2 when 
the supergravity Kahler potential takes a particular form, e.g. the no-scale form flT5| . 

Having defined our supersymmetric axion model, we discuss the full 4x4 axino-neutrino 
mass matrix: ^ 

-m AB v A v B (10) 

where A,B = s, e, //, r and v s = a with m ss = rria. We will work in the field basis in which 
li' i LiH 2 (i = e, fA, t) in W e R are rotated away by an appropriate unitary rotation of H 1 and 
Li. It is straightforward to see that the axino-neutrino mixing mass is given by 

/pg UO-V V 600 GeV/ \ r ! 



where q = /4//io f° r K i n Weff before fi' i L i H 2 are rotated away. Note that this axino-neutrino 
mixing survives under the unitary rotation eliminating fi' i L i H 2 since Lj and ifi have different 
£/(l)pQ-charges. This charge difference is also responsible for the suppression of _R-parity 
violating couplings as well as the weak-scale value of iiq. The 3x3 mass matrix of active 



neutrinos is induced by -R-parity violating parameters [ITHj. At tree-level 



a 



where M a denote the gaugino masses. The sneutrino VEV's (z/j) are determined by the bi- 
linear i?-parity violations in the SUSY-breaking scalar potential: m 2 LiHi LiH\ + B[LiH2- In 
our model, nonzero values of rn 2 L H and B[ at the weak scale arise through renormalization 
group evolution (RGE), mainly by the coupling X' i33 yb where y^ is the 6-quark Yukawa cou- 
pling ||17||. Moreover, BHiH 2 arises also through RGE which predicts a large tan/3 ^40 — 60 



We then find 



mij „ 10 -V M m) eV (13) 
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where t = \n(Mx/mj) / ln(10 3 ) for the slepton mass mj. Here we have taken mj pa 300 
GeV and //q ~ 2mj which has been suggested to be the best parameter range for correct 



electroweak symmetry breaking [18|] . There are also bunch of loop corrections to m^ from 
.R-parity violating couplings ||19|| , however in our case they are too small to be relevant. 

Let us now see how nicely all the neutrino masses and mixing parameters are fitted in 
our framework. The analysis of Ref. J?| leads to the four parameter regions, R1-R4 of Table 
1, accommodating the LSND with short baseline results. In our model, Eqs. @ and ( |TT|) 
can easily produce these LSND values of the largest mass eigenvalue m 4 pa m ss = rria ~ 1 
eV and the mixing matrix elements ~ rn is /m ss ~ 0.1. The effective 3x3 mass matrix 
of active neutrinos is then given by 



eff 

™>ij = m ij 



m is m js /m ss . (14) 



Upon ignoring the small loop corrections, this mass matrix has rank two, and can be written 

as 

mtf = m x XiXj + myjijj (15) 

where %, and f j are the urr h vectors in the direction of m ls and (i/j), respectively. Remark- 
ably, the mass scale m x ~ (mi S /m ss ) 2 m ss ~ 10~ 2 eV gives the right range of the atmospheric 
neutrino mass. Eq. ( |T3"D shows that m y is also in the range of 10 -2 eV, so m eff would be 
able to provide the right solar neutrino mass unless Am^ <C 10~ 4 eV 2 . Note from Eq. (|j) 
that the typical size of e«, X ijk , X' ijk is around 10~ 6 for f PQ pa 10 10 GeV and M* pa 10 16 GeV. 
To discuss the consequeces of the mass matrix (|i~5|) on the atmospheric and solar neutrino 
oscillations of three active neutrinos, we diagonalize it to find the mass eigenvalues 

^2,3 = - (m x + my±^J (m x + m y cos 2 2£) 2 + m 2 sin 2 (16) 

and the 3x3 mixing matrix 

U = (z T , w T c e - x T s e , ~w T sg + x T c e ) . (17) 

where z = x x y/|x x y|, w = x x z/|x x z|, and tan26 l = m y sin 2£/(m x + m y cos 2£) for 
cos^ = x-y. The Super-Kamiokande data [^] combined with the CHOOZ result [^TJ imply 
that t/ 2 3 pa U^ 3 pa 1/2 and ?7 2 3 -C 1. The solutions to the solar neutrino problem can have 
either a large mixing angle (LA): ~ £/ 2 2 ~ 1/2, or a small mixing angle (SA): pa 1. 
This specify the first column z T of U as z pa (l/\/2, — 1/2, 1/2) for LA and z pa (1,0,0) 
for SA up to sign ambiguities. Since x • z = 0, the pattern z pa (1,0,0) implies x e pa 0. 
This leads to a too small ?7 e 4 pa m es /m ss < 10 -2 , so the SA solution is not allowed within 
our model. Among various LA solutions to the solar neutrino problem, only the large-angle 
MSW solution with Am 2 oi ~ 10 _4 eV 2 can be naturally fitted since m x pa m y ~ 10 -2 eV in 
our scheme. It is remarkable that J'pq ~ 10 10 GeV and pa Mqut lead to the right size of 
.R-parity violation yielding the desired values of m is and my also for the atmospheric and 
solar neutrino masses. 

To see the feasibility of our whole scheme, we scanned our parameter space to reproduce 
the allowed LSND islands R1-R4 of Table 1 together with the following range of atmospheric 
and solar neutrino parameters: Am| x = (1 — 8) x 10~ 3 eV 2 , Am 2 x = (0.1 — 8) x 10 -4 eV 2 and 
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Am^eV 2 ) 


\U ei \ 


IM 


Rl 


0.21-0.28 


0.077-0.1 


0.56-0.74 


R2 


0.88-1.1 


0.11-0.13 


0.15-0.2 


R3 


1.5-2.1 


0.11-0.16 


0.09-0.14 


R4 


5.5-7.3 


0.13-0.16 


0.12-0.16 



TABLE I: Allowed regions for the LSND oscillation. 
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FIG. 1: Scatter plots on (a)m y vs. m x , (b)cos£ vs. 
parameters for R2. 
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x reproducing the correct oscillation 



tan 2 #23 = 0.33 — 3.8, tan 2 9yi = 0.2 — 3.0, tan 2 9i 3 < 0.055 [|TJ. Our parameter space consists 
of m ss ,m is , \' i33 yb whose values are centered around 1 eV, 0.1 eV, 10~ 6 , respectively. For Rl 
and R4, we could find some limited parameter space, however they need a strong alignment 
between x and j and also a large cancellation between m x and m y . On the other hand, 
R2 and R3 do not require any severe fine tuning of parameters, although the bi-maximal 
mixing is obtained by some accident. We provide in Fig. 1 the scatter plots on the planes 
of (m y ,m x ) and (cos£,m x ) for the LSND island R2. The scatter plots for R3 have similar 
shape. The full details of the analysis will be presented elsewhere JT2] |. 

To conclude, we have shown that the 3+1 scheme of four-neutrino oscillation can be 
nicely obtained in supersymmetric model endowed with the PQ solution to the strong CP 
problem and supersymmetry breaking mediated by gauge interactions. 
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